Flux-noise spectra around the Kosterlitz-Thouless transition 
for two-dimensional superconductors 
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The flux-noise spectra around the Kosterlitz-Thouless transition are obtained from simulations of 
the two-dimensional resistively shunted junction model. In particular the dependence on the distance 
d between the pick-up coil and the sample is investigated. The typical experimental situation 
corresponds to the large-d limit and a simple relation valid in this limit between the complex 
impedance and the noise spectra is clarified. Features, which distinguish between the large- and 
small-d limit, are identified and the possibility of observing these features in experiments is discussed. 
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I. INTRODUCTION 



Spontaneously created vortices drive the transition be- 
tween the superconducting and normal state for thin-film 
superconductors and twa-dimensional (2D) Josephson- 
Junction arrays (JJA's)IH This means that the physics 
of the vortices is responsible for the characteristic fea- 
tures in a region around the transition. One manifesta- 
tion of this is the static characteristics of the phase tr. 
sition which is of the Kosterlitz-Thouless (KT) type 
Another manifestation is the dynamical features around 
the transition. These are reflected in the flux-noise spec- 
tra and the complex impedance. In the present paper 
we investigate the flux-noise spectra through computer 
simulations of the resistively shunted Josephson junction 
(RSJ) model on a 2D square lattice. In particular we 
clarify the relation between the flux-noise spectra and 
the complex impedance. 

There have been a number, aLrecent experimenta 
as well as theoretical studiesQlij dealing with the flux- 
noise spectra. The typical experimental setup measures 
the fluctuation of the magnetic flux through a -pick- up 
coil situated at a distance above the sample. ol3 Many 
simulation studies on the other hand have measured the 
vorticity fluctuation associated with a fixed area of the 
sample itself.ETEj It has been assumed that this would 
roughly correspond to the magnetic flux spectra of the 
experiments. However, in the present investigation we 
show that there are significant differences. The typical 
experimental situation corresponds to the limit of large 
distance between the pick-up coil and the sample. In 
this limit there exists a simple relation^ between the flux- 
noise spectra and the complex impedance of the sample, 
which we here verify both directly from the simulations 
and through analytic calculations. 

In Sec. O we describe how the flux-noise spectrum is 
obtained from simulations of the RSJ model. Section III 
clarifies the relation between the flux-noise spectrum and 
the complex impedance. The results from the simulations 
are described and discussed in Sec. IV. Particular atten- 



ments of flux noise and the complex impedance. Finally 
Sec. M contains some concluding remarks. 



II. RSJ MODEL AND FLUX NOISE 
A. RSJ model and numerical method 

In our simulations we use the 2D RSJ model on a 
square lattice with periodic boundary conditions (PBC). 
This is usually assumed to be a good model of a Joseph- 
son junction array (JJA). In particular, from the point 
of view of vortex fluctuations and vortex dynamics, it is 
expected to have the same physics as a thin supercon- 
ducting film, as well as a JJA.EJ There might however be 
some differences, in the level of vortex fluctuations around 
the transition.113 

The RSJ model incorporates the condition of the local 
current conservation and the equations of motion can be 
written as:u 



= - Yl Gi 3 Yl ( sin ^ k + ^ k ) ' 



(1) 



where Gij is the square lattice Green function, the primed 
summation is over the four nearest neighbors of the site 
j, and <j)jk = Oj — 9k with the phase 9j of the complex 
order parameter at site j. Here we measure time t in 
units of h/2eRI c , where R is the shunt resistance and I c 
is the critical current of a single junction. The thermal 
noise current r\jk in units of I c satisfies the conditions 
(Vij(t)) = and 



(Vij(t)Vki(0)) - 2T(5 lk 5 jl - 6 a 6 jk )6(t), 



(2) 



tion is given to the implication for experimental measure- 



where (• • •} is the ensemble average, and the tempera- 
ture T is in units of the Josephson coupling strength 
J = HI c /2e. 

The RSJ model may be used to calculate the cur- 
rent distribution in the limit of a large perpendicular 
penetration length A = &oc/4ir 2 I c , where $o is the 
flux quantum^ This is typically the case for a 2D su- 
perconductor.S In this limit one may replace the gauge 



invariant phase difference 4>jk = Oj — Ok — Ajk, where 
Ajk = |r J- A(r) • d\ is the line integral of the vector 
potential, with cjjjk = Oj — Ok because the coupling to the 
electromagnetic self-field is in this limit so weak that it 
has little influence on the fluctuations of the supercur- 
rent (we are here considering the situation without an 
external magnetic field) £l 

We measure the flux-noise spectrum S(oj) defined by 
the Fourier transformation: 



from the simulation data. However, when presenting the 
data we have for convenience used an extrapolation to 
large t based on an expected large-i behavior (this ex- 
trapolation does not change the behavior in frequency 
range u> > 0.016). From the ergodicity of the system, we 
can change the ensemble averages of the form (O(t)O(O)) 
in Eqs. @ and (0) to the average (0(t + t')0{t')) t , over 
time if , and averages over more than 10 7 time steps were 
typically performed. 



S(w) 



dte iut S(t), 



(3) 



where S(t) is the time-correlation function of the mag- 
netic flux $(£) through a pick-up coil: 



S(t) = <*(t)*(0)> 



(4) 



(see Sec. [IB for the calculation of $). In addition, we 



calculate the dynamic dielectric function \/dxj) of the 
vortices in the Coulomb-gas analogy given bycj 



Re 



Illl 



1 



e(co) 
1 



1 



2ttujT cg f°° 



T 2 



dt Bin wtG(t), (5) 



c( W ) 



2ttluT cg 
T 2 



dt cos ujtG(i), 



(6) 



where T CG = T/(2ttJ(cos4>)) is the Coulomb-gas tem- 
perature, and the time-correlation function G(t) is de- 
fined by 



G(t) 



1 



L 2 
F(t) = J2 sin 



(F(t)F(0)), 



'<) ■ 



(7) 
(8) 



(iJU 



and the sum is over all links in ^-direction. 

The dynamic dielectric constant l/e(w) is related to 
the-conductivity <j{uj) and the complex impedance Z(u>) 

b y y 



a(u) = Z -1 (w) = 1 



T 



1 



i2irojT CG e(w) ' 



(9) 



In the numerical simulation of an L x L array (in most 
cases L = 64 and occasionally L = 128 are used in the 
present paper), we use the periodic boundary condition 
for the phase variables, i.e., #i+Lx = #i+Ly = #i, and 
the thermal noise currents are generated from the uni- 
form probability distribution. For the time integration of 
equations of motion in Eq. p), we use the Euler method 
with the discrete time step At = 0.05. In practice we 
have calculated S(t) up to a i max beyond which S(t) be- 
came so small that the simulations could not be con- 
verged well enough to obtain further information. In our 
simulations this turned out to be i max w 100 for T > 1.10 
and i max s=s 400 for lower temperatures. This means that 
we could not reach frequencies below lj < 0.016 directly 



B. Flux noise and the dipole approximation 

As mentioned above, previous numerical simulations&EI 
usually calculate the noise spectra from the fluctuation of 
the vorticity defined by the directional sum of the gauge- 
invariant phase difference around each plaquette. The 
fluctuations of the total vorticity over an area of the sam- 
ple is then used to estimate the flux-noise spectrum. On 
the other hand in experimentstTcl the fluctuations of the 
magnetic flux penetrating a pick-up coil situated at a 
distance above the sample is measured. The relation be- 
tween the vorticity noise and the flux noise has so far not 
been studied in detail and is the subject of the present 
paper. 

We consider three distinct cases: The first one is the 
fluctuation of the magnetic flux associated with the vor- 
tices on a fixed area A of the sample. The vortices dfe 
scribe the rotation of the supercurrent on the sample.EJ 
This means that the magnetic flux associated with an el- 
ementary plaquette is proportional to the rotation of the 
supercurrent around the plaquette. For the RSJ,-mQd<;l 
the magnetic flux for a plaquette is then given bj 



IK 



n(r) = 



T 



C'G 



T 






"y > 



(10) 



E 



in units of the flux quantum 0o -ell Here r is the central 
position of the plaquette and the summation is around 
the plaquette p. Consequently the total magnetic flux as- 
sociated with the rotation of the supercurrent at a given 
time t is 



$(i) 



d 2 r n(r,t), 



where the integral over position denotes the sum over all 
elementary squares inside the area A. 

The second case is the fluctuation of the vorticity as- 
sociated with a fixed area A of the sample. The vorticity 
of an elementary plaquette is given bycl 



v(r) 



2^' 



(11) 



where the phase difference </>y is restricted to the interval 
— 7r < <f>ij < tt. The total vorticity V(t) associated with 
the area A is consequently 



V(t) = / d 2 rv(r,t). 

J A 

The third case corresponds to the experimental situ- 
ation where one measures the magnetic flux through a 
pick-up coil situated a distance d from the sample. Note 
that the first case corresponds to the case when d = 0. 
However, in typical experiments d is,ajnacroscopic length 
(typically of the order of 0.1 mm).cTEl 

In this section, we focus on the third case and use 
a dipole approximation to obtain an expression for the 
magnetic flux through a pick-up coil at a distance d from 
the sample where d> A. At this distance the magnetic 
field distribution from a vortex is of dipole form. In the 



continuum limit only the circulation of the supercurrent 
around a closed loop which encloses a vortex core region 
gives a finite value. This value can be positive or negative 
but has the same magnitude for all closed loops which 
enclose vortex cores E This means that one can estimate 
the magnetic field at d 3> A from the circulation of the 
supercurrent around small closed loops which cover the 
area of the superconductor by associating each of these 
loops with a dipole field where the strength of the dipole 
moment is proportional to the circulation of the super- 
current. This approximation is readily carried over to 
the array. The circulation of the supercurrent around a 
plaquette in units of I c is given by [see Fig. |I](a)]Ea 
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(12) 



where </>(r',r") = 8 r > — 9 T n and 6 r > is the phase of the complex order parameter at site r'. The contribution to 
the magnetic field at large distances from this circulation can then approximately be described in terms of a dipole 
moment m r = m r z proportional to the circulation. The magnetic field at the observation point x due to this vortex 
dipole moment m r is given byt3 



Bfxl 



E 



3e rx (e rx • m r ) m r 



(13) 



where e rx = (x — r)/ |x — r | is the unit vector in the direction of x — r and the summation is over all dual lattice points 
(central positions of plaquettes). The magnetic flux through the pick-up coil [the shaded area in Fig. 0(b)] separated 
by distance d from the array is given by the surface integral: 



$= / B(x) -da. 

./coil 



(14) 



We consider the case d > a so that the magnetic field inside the pick-up coil does not vary much on a microscopic 
length scale a. We may then replace the integral Eq. dlj) by the discrete summation: 



$ « Y^ B(r' + dz)-z 

r'Elxl 

•sr-^ sr-^ 3n rr '(n rr ' • m r ) — m r 



r'eixZ reLxL 



|r' + dz — r|- 



v^ v^ 3d 2 - |r' + dz- if 

E E — — - ^ TO * 



r'Elxl reLxL 



|r' + dz — rl 5 



(15) 



where x is decomposed into x = r' + dz, and the two 2D vectors r and r' denote positions on the array. The summations 
^2 T> and ^ r are performed on the I x I pick-up coil and the LxL whole array, respectively [see Fig. Mb)] . Since Eq. ( |15| ) 
contains 0(^ 4 ) terms (in this work we choose I = L/2), the calculation of the magnetic flux in this way requires most 
of the computer time in the numerical simulations. This time-consuming part is avoided in the approximate scheme 
we use in our simulation to obtain the results described in Sec. 



IV 



We use the following approximate scheme: The whole array is divided into stripes formed by elementary plaquettes 
which enclose the midpoint of the pick-up coil as the sides of a square. This is illustrated in Fig. 0(c). Each such 
collection of elementary plaquettes are denoted by S n where 2n— 1 is the number of plaquettes of each of the four stripes 
forming the sides of the square. The magnetic flux in Eq. ( jig ) is expressed as the summation of the contributions 
from each S n : 



$ = J2 1nM n 



(16) 



where M n is the summation of the vortex dipole moments for the plaquettes forming S n and q n is the appropriate 
weight factor: 



?™ = ( E TOr E 



Vr£S„ r'eZxi 



3d 2 - |r' + dz-r\ 
|r' + dz — rl 5 



M re = ^ 



r£S„ 




1 



(17) 
(18) 



We now assume that m r in Eq. ([l7|) can to good ap- 
proximation be replaced by the average value for each 
S n , i.e., m r = J2res rn r/ j ^-n where A n is the number of 



plaquettes contained in S n . This simplifies Eq. (17) to 



Qn 



1 



E E 

r6S„ r'GZxZ 



3d 2 



rfz 



(19) 



The magnetic flux in Eq. (jig ) can by aid of Eqs. ( (L8[) and 
(119) be computed in O(Lr) operations, since q„ within 
this approximation is a purely geometric quantity which 
is independent of time. In Fig. H we compare the flux 
noise S(t,d) = ($(i, d)$(0, d)) from the full calculation 
in Eq. ( jig ) and the ^-approximation in Eqs. ( |l6| ) and 
( |l9| ) for a 32 x 32 array with a 16 x 16 coil size and the dis- 
tance d = 5. It is clearly shown that the approximation 
made in Eq. ( J19J ) is indeed a very good approximation. 

Figure shows q n as a function of the linear size 
(2n - 1) of ,S„ for d = 0.1, 10, and 20 (a 32 x 32 pick-up 
coil is used for a 64 x 64 array). For very small values 
of d, q n becomes a step function where only S n 's inside 
the pick-up coil contribute to the magnetic flux. On the 
other hand, as d is increased, it is clearly seen that there 
is a significant contribution to the magnetic field caused 
by the S n 'a outside the pick-up coil area. In previous 
numerical studies of the flux-noise spectra the magnetic 
flux has usually been approximated by the vorticity in- 
side the pick- up coil area (the second case mentioned in 
the beginning of this subsection) nu This approximation 
hence does not take the contributions from S n 's outside 
the pick-up coil area into account and in this sense it cor- 
responds to d = 0. One conclusion from the present work 
is that for a more precise comparison with experiments 
one should instead consider the opposite case of large d. 



The results of our simulations are presented in Sec. IV 



In the following section we elucidate the relation between 
the flux-noise spectrum and the complex impedance, or 
equivalently the complex conductivity. 



III. FLUX NOISE AND CONDUCTIVITY 

In this section we consider for simplicity a 2D su- 
perconductor in the continuum limit so that, compared 



to the previous section, the limit a — » is implied in- 
stead of a = 1. The magnetic flux associated with the 
area d 2 r around r is then n(r)d 2 r [see Eq. (O)]. In 
the Couipmtt-gas analogy of vortices, n(r) is the charge 
densityElo'EEl The charge density correlation function is 
given by c(r,t) = (n(r, t)n(0, 0)). We will first relate 
c(r,t) to the dielectric function l/e(k, u>) and the con- 
ductivity a(u/) of the superconductor: The charge den- 
sity correlation function c(r, t) is related to the charge 
density response function g(r, t) by 



Im[g(k,uj)] = 



2T CG 



c(k,w), 



(20) 



where g and c denote the Fourier transforms. The di- 
electric function |J/e(k, cj) is given by the usual linear 
response relation! 



e(k, u) 






(21) 



We define l/e(oj) = l/e(k = 0,w) which means that 
Eqs. ( p0[ ) and (Ell) for k = corresponds to Eqs. (|) and 
(|) for the RSJ model. From Eqs. (§), ©, and @ we 
obtain a relation between the charge density correlations 
c(r,t) and the real part of the conductivity &(lo): 



and 



Re[<rM 



Im 



T 



2ttT cg l 



-Im 



1 



e(w) 



e(w) 



ituj c(k,w) 

'tcgJ^ — p— 



(22) 



(23) 



Next we relate the charge density correlation function 
c(r, t) to the flux-noise spectrum. From Eq. (15) we have 
that the magnetic flux measured by the pick-up coil is 



$ = I B z (r)d 2 r, 

J coil 

where the integral is over the area covered by the coil. 
The magnetic field B z (r) can be expressed as 

B z (r) = J f(\r'-v\,d)n(r)d 2 r', 



the r'-integration is over the whole 2D plane, and from 
Eqs. © and (||), we have 



rpCG 



(r 2 + d 2 )t 



This means that the flux-noise spectrum 5(t) 
($(t)<f>(0)) is given by 



S(t) = / d 2 r 



coil 



coil 



dV / d 2 r" / dV'7(|r"-rU)(n(r,t)n(r / ,0))/(|r'-r'"|,d) 



(25) 



We can now use the convolution theorem and express S(t) in terms of the Fourier transforms of /(r, d) and c{r, t) 
(n(r,t)n(0,0)>, i.e., 



sw =l/ r L, dV /(^ e< 



ik-(r-r'; 



|/(M)| 2 c(M). 



Taking the Fourier transform of S'(i) and using the con- 
nection between c(k,u>) and l/e(k, w) given by Eqs. (|20| ) 
and (|2l|) yields 



S(w) = - / dfcF(fc, d)Im 

Jo 



f 



e(k,ui) 



(26) 



where 



F(k,d) 



(2tt)^w 7 coil J coil 



(27) 



and /(fc, d) is the Fourier transform of /(r, <i) which de- 
scribes the spreading of the magnetic field. Within the 
dipole approximation of Eq. (|2J) f(k, d) is given by 



f(k,d) 



T 

rpCG 



2itke 



-kd 



(28) 



The extreme case d — is outside the dipole approxima- 
tion and is given by 



/(*) 



T 

rpCG' 



(29) 



The important point to note is that the flux-noise spec- 
trum S(lu) is directly related to the response function 
lm[l/e(fc,a;)] through a function F(k,d) which contains 
all the information of the spreading of the magnetic field 
above the superconductor as well as the geometry and 
position of the pick-up coil. 

In Ref. it was suggested that the relation between 
the flux noise S(uj) and Im[l/e(A;, w)] could be further 
simplified to 



S(lu) oc — 



Im 



Le(0,w) 



or equivalcntly 



S(u>) «Re[(T(w)]. 



(30) 



We will here show that for the typically experimental sit- 
uation this proportionality between the conductivity and 
the flux noise is indeed valid. 

In order to establish this we assume for simplicity that 
the pick-up coil is circular with radius R. In this case it 
is possible to obtain an explicit expression for -F(fc, d) in 
Eq. (§7|), i.e., 



F(k,d) 



2T cg 



i? 2 /fc|/(A:,d)! 2 [J 1 (fci?)] 2 



(31) 



where J\ is the Bessel function of order one. Let us 
first consider the limit in which the dipole approximation 
Eq. (|2J) is valid. This limit implies that the distance d to 
the pick-up coil is sufficiently large compared to the rele- 
vant microscopic lengths. For a Josephson junction array 
the microscopic length is the lattice constant a which is 
typically of the order 1-10 /im whereas for a continuum 
superconductor it is the size of a vortex core typically 
given by the Ginzburg-Landau coherence length £ and 
is of the order of 100-1000 A close to the transition. In 
addition A has to be much larger than a for a Joseph- 
son junction array and £ for a,s«perconducting film. The 
typical size of d is 100 /im.ETEl So in practice the vor- 
tex dipole approximation may be expected to be valid 
for superconducting films but the validity for Josephson 
junction arrays may be more questionable. When the 
vortex dipole approximation is valid the noise spectrum 
is given by [compare Eqs. (26), (E8), and (31)] 



Sr(oj) 



lir 2 T 2 R 2 
TCGlu 



dkk 6 e-^ d [Ji(kR)Yhn 



i{k,uj) 



(32) 



One notes that k- values much larg er than 1/d will not 
contribute to the integral in Eq. (|3^) because of the factor 
e~ 2kd . This means that if d is sufficiently large compared 
to the relevant microscopic length, then Im[l/e(fc,ct>)] 
can be replaced by lm[l/e(0,<x>)], demonstrating that 
under these conditions S(u) is indeed proportional to 
lm[l/o;e(0,a;)]. In the present case of a circular pick-up 
coil we explicitly find 

C, 



Sr(u>) = Im 



e(0,o;) 



where the proportionality constant C is given by (after 
changing the integration variable to x = kd) 



C 



8ir 2 T 2 R 2 
~T Ug ~1F 



dxx e 



■2.1- 



[Ji{xR/d)\< 



In the limit R/d <C 1 this reduces to 



C 



T 2 R 4 157T 2 



and in the limit R/d ^ 1 to 



C 



T 2 2ttR 
T CG d 3 " 



(33) 



(34) 



We conclude from Eqs. ( p3[ ) and (pil) that the propor- 
tionality constant C will always contain a temperature- 
dependent factor T 2 /T CG and a dependence on the size 
of the coil. This coil-size dependence reflects the relative 
magnitudes between the coil size and the distance from 
the sample: When R is much larger than d, the noise am- 
plitude is proportional to the perimeter of the coil 2irR, 
and, when it is much smaller, it is proportional to square 
of the coil area i? 4 . n Ia typical experiments R is usually 
much larger than daU 

In some previous simulations of the flux-noise spec- 
trum, based on the XY models, one has approximated 
the flux-noise spectrum from the fluctuation of the total 
vorticity for a finite area of the model.™ This implies two 
differences in relation to the above dipole approximation: 
First of all it corresponds to d = and consequently to 
a constant / [see Eq. (p9|)]. Secondly it corresponds to 
changing the magnetic flux defined by Eq. (HI) to vor- 
ticity defined by Eq. (|ll]). Let us first consider the first 
change by itself: The case of a circular area with radius 
R then corresponds to the flux noise [compare Eqs. (26), 
©, and @] 



Sr(uj) = 



2T 2 

rpCG J 



OJ 



R 2 / dfcfc[Ji(fci?)] 2 Im 



1 



_i(k,uj) 



(35) 



This means that, in this case, the flux-noise spectrum 
depends on all the /c-valucs of Im[l/e(fc,u;)] and is not 
proportional to lm[l/e(0,cj)]. For example the leading 
large- R dependence of Eq. (pS) is 



Sr(uj) cx / dklm 

v Jo 



e(k,uj) 



So in this limit the flux noise is proportional to the 
perimeter of the pick-up area just as for the large-d case, 
but instead of singling out the k — O-contribution all 
/c-values contribute. The proportionality between the 
flux noise S(tu) and lm[l/u;e(0, lu)] cx Re[<r(aj)] can be 
tested by experiments since both S(to) andrihe conduc- 
tivity a{uo) can be independently measured.tl 

The change from magnetic flux [defined by Eq. jiq)] 
to vorticity [defined by Eq. (Ill])] influences the flux-noise 
spectrum in an additional significant way: Now the cross- 
ing of a vortex over the perimeter of the pick-up area is 
described as a discrete ±27r-change of the total vorticity 
of the pick-up area. The corresponding spectrum hence 
corresponds to a random walk of discrete events over a 
sharpj-boundary This results in a w~2-tail of the spec- 
trumEl However, this condition does not correspond to 
the experimental situation where the pick-up coil does 
not have a sharp boundary, is at a distance d from the 
sample, and, most importantly, the magnetic field from 
a vortex is spread out. 

In the next section we present numerical results of the 
flux-noise spectrum and its relation to the conductivity. 



IV. SIMULATION RESULTS AND 
EXPERIMENTAL IMPLICATIONS 



A. Comparison with previous works 



We first relate oux -simulations results to earlier ones 
for the RSJ model.Em These earlier simulations calcu- 
lated the noise spectrum of the vorticity [see Equity])] 
over a fixed area of the systems (the d = 0-case) ImAs 
explained in the previous section, this corresponds to dis- 
crete-events over a sharp boundary and implying a u>~ 3 ' 2 - 
tailEJ Such a tail has indeed beea found in Ref. and is 
also verified in our simulations HH This is apparent from 
Fig. [| which displays our data for a 64 x 64 array with a 
pick-up area of size 32 x 32 at T = 1.1; the lower data- 
set shows the vorticity-noise spectrum and the slope is 
—3/2. However, if we instead calculate the magnetic- 
flux- noise spectrum [see Eq. (JL0|)], then the crossing of 
magnetic flux over the perimeter is not a discrete event. 
This means that there is no obvious reason for a w~ 3 / 2 - 
tail and nor do we find any such tail in the simulations, 
as is also apparent from Fig. [|. The exponent in Fig. [| 
for this case is instead close to —2 as seen from the upper 
data-set in Fig. ||.c3 



B. Flux- noise spectra with a finite distance between 
pick-up coil and array 



As mentioned in Sec. [II the typical experimental setup 
measures the magnetic- flux-noise spectrum with|-a.-finite 
distance d between the array and the pick-up coiLuEl Fur- 
thermore the typical experimental setup corresponds to 
the large-d limit where the noise spectrum S(u)) is pro- 
portional to the real part of the conductivity (t(uj) (see 
Sec. 0). 

In our simulations we investigate the flux-noise spec- 
tra as a function of the distance d to the pick-up coil. 
Figure 3(a) shows the flux-noise spectra calculated from 
Eqs. (|), (4|) with the mag netic flux given by Eqs. (16), 
(|l8|), and (^j) (see Sec. IIB for details). The data sets are 
shown as uoS{bj, d) against w in a log-log plot. The verti- 
cal scale is adjusted in order to compare the shapes of the 
curves. One notices that the spectra for the different d's 
all approach w _1 for large u> and can be collapsed to a sin- 
gle curve in this large-w region by a vertical adjustment. 
For small to the curves becomes linear with WnWrhich re- 
flects a constant part (white noise) of S(uj, <i).ErO£m As d 
increases the peak of the ujS(oj, d) curves moves to the left 
and the peak height increases. The uppermost curve in 
Fig. ||(a) is |Im[l/e(w)]| (full curve in the figure) and, as d 
is increased, the flux-noise spectrum ujS(uj, d) approaches 
this uppermost curve. This verifies that for large d one 
has the simple connection S(uj) oc |Im[l/we(a;)]| as dis- 
cussed in Sec. [Hl| 

Figure 0(b) shows that the characteristic frequency 
given by the peak position in Fig. 0(a) decreases with 
increasing d. In the limit of large d the characteristic 
frequency of S(uj,d) agrees with the characteristic fre- 
quency of l/e(w). Thus in this limit both the shape and 
the characteristic frequency of S(u,d) and Im[l/u;e(w)] 
are the same. This proportionality between the flux-noise 
spectrum and complex conductivity can be tested exper- 
imentally and indeed seems to be borne out.cl 

The fact that the flux-noise spectrum in the large- 
d limit is proportional to real part of the conductivity 
means that the characteristic features of the conductiv- 
ity are reflected in the flux-noise spectrum. In case of 
a 2D superconductor the dynamical features of the con- 
ductivity <j(lu)— pc — l/iuje(uj), are well described by the 
response formtil 
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which catches the dynamics of vortex fluctuations in a re- 
gion around the KT transition. From Eq. ( fi7\ ) one notices 
that the peak of |Im[l/e(w)]| occurs at the characteristic 
frequency loq and that the peak height is 1/ire. Above 



the KT transition 1/e increases only weakly with in- 
creasing temperaturej-apd approaches unity for somewhat 
higher temperatures Jiil For the flux-noise spectrum this 
means that S{ui ) ex T 2 /T CG iuj . Now T/T CG ex p (T) 
where po(T) is the bare superfluid density which de- 
creases slightly with temperaturetl whereas T increases 
so that also T 2 /T CG depends only weakly on T. This 
means that to good approximation the flux-noise spec- 
trum for different temperatures above the KT transition 
should have a common tangentoc 1/w which goes through 
all the points S(ujo(T), T). This means that the common 
tangent in a log-log plot has the slope —1. This feature 
is illustrated in Fig. @ which demonstrates the existence 
of a common tangent with the slope close to — 1 both di- 
rectly for S(ui,d) and for Re[<r(a;)] oc —lm[l/uie(w)]. The 
existence of a common tangent for the flux-noise spectra 
with the slope —1 can be readihj— tested in experiments 
and seems to be well borne out.BE3 

One should, however, notice that the argument for a 
common tangent with slope —1 does not single out the 
response form given by Eqs. ( pq ) and (p7|). In fact the 
reasoning is also valid in a region above the KT transition 
where the response of the vortex fluctuations is given by 
the conventional Drude response form 



and 



Re 



Im 



eH 



(«) 



ewHffi 



1 LOLUq 



(38) 



(39) 



which gives the peak height l/2e. One expects that 
the response form Eqs. ( |36| ) and ( p7[ ) describes the re- 
sponse from the vortex pairs in a region just above the 
KT transition whereas the conventional Drude response 
is obtained for higher temper.atures where the response 
is dominated by free vortices £3 How wide these regions 
are depend on the details: for a real thin superconduc- 
tor the,-wortex-pair response seems to dominate in a wide 
region. til However, for the 2D RSJ model on a square 
lattice, which we are using here, the vortex-pair domi- 
nated region above the KT transition is narrow and the 
Drude response dominates in a broader region above the 
KT transition. Thus the data shown in Fig. are pre- 
dominantly Drude-like. A practical way of determining 
which response type is at hand is to measure the complex 
impedance and determine the peak ratio [i.e., the ratio 
Re(<r)/Im(cr) at the peak position ujq of Re(a/ui)]; for the 
vortex pair dominated response, Eqs. ( pq ) and (|37J), this 
ratio is 2/tt w 0.64 and for free vortex Drude response, 
Eqs. (||) and ©, it is unityN 

The essential point here is that, because the flux-noise 
spectrum for a large d is proportional to Re[cr(o;)], the 
noise spectra at a sequence of temperatures above the KT 
transition should in a log-log plot have a common tangent 
with slope —1. We can substantiate this claim further by 
simulating the noise spectra for a small d where the flux- 
noise spectrum is not proportional to Re[er(aj)]. In this 



case there is no particular reason for a common tangent 
with any slope and, as apparent from the simulation re- 
sults in Fig. 7], no such common tangent can be fitted to 
the data. 

One may also notice from Fig. that both S(ui,d) 
and Re[cr(w)] have intermediate regions with oj -1 - 5 fol- 
lowed by a w~ 2 -tail for even larger ui. However, such an 
intermediate-^ -1,5 region appears to be less discernible 
for the sma\\-d case when S(oj, d) is not proportional to 
Re[er(w)], as is apparent from Fig. R. 



C. Flux- noise spectrum below KT transition 



tion is at T w 0.9 and already at T = 0.85 the amplitude 
of u;S(uj,d) has dropped dramatically compared to the 
almost constant amplitude above the KT transition, as 
is apparent from Fig. |S|(b). 

Finally, there is in Fig. [81(b) an indication that the 
curves develop a plateau as the KT transition is ap- 
proached from above (compare the curve for T — 0.95). 
Such a plateau would suggest that S(u>,d) is propor- 
tional to 1/uj in an intermediate region just above the 
KT transition. The same development of a plateau can 
be anticipated in Fig. |](a) for |Im[l/e(o;)]| and has also 
been found for the XY model with thc.-time-dependent 
Ginzburg-Landau (TDGL) dynamicsHEl 



Next we investigate what happens as the temperature 
is decreased towards the KT transition and below. Fig- 
ure g shows data for Im[l/e(w)] and ujS(lo, d = 20) over 
a wider range of temperatures (the data for T > 1.1 
arc the same as in Fig. g]). Again one observes that 
both quantities behave in precisely the same way over the 
whole temperature range, verifying that they are indeed 
proportional to very good approximation. Next one ob- 
serves that the characteristic frequency luq (the frequency 
of the peak position) decreases as the KT transition is 
approached from above. This suggests a .critical slowing 
down at the KT transition to ujo = O.ErE-3 As the tem- 
perature passes through the transition the characteristic 
frequency starts to increase agap,(the full curves in Fig. || 
are below the KT transition) Bt3 

The argument for a common l/w-tangent for the flux- 
noise spectra above the KT transition is related to the 
fact that the peak height for |Im l/e(u)] is 1/ire (l/2e) for 
the vortex pair response, Eqs. ( 36[) a nd (|37| ) [free vortex 
Drude response, Eqs. 



) and (p9|)] together with the 
fact that 1/e only increases very weakly with T above 
the KT transition and approaches unity for somewhat 
higher T . Figure [S| shows this weak increase in a re- 
gion above the KT transition; the Drude value 1/2 is 
approached roughly like w which explains the discrep- 
ancy between exponent 1 and the value 0.9 found in 
Fig. U Thus the existence of a common tangent 1/uj 
hinges on the weakness of the temperature dependent 
factor T 2 /T CG (T)e(T) which in turn depends somewhat 
on the details of the system. However, since the temper- 
ature dependence of loq is dramatic just above the KT 
transition, the common tangent should to good approx- 
imation exist at least in a limited region above the KT 
transition. 

Below the KT transition there are no free vortices and 
the response is given by the vortex pairs Eqs. (pq) and 
(p7|) .E-3 However, in this case the factor 1/e w 1 — l/e(w = 
0) decreases rapidly towards zero as-the temperature is 
decreased below the KT transitiono This means that 
while the characteristic frequency rapidly increases, as 
the temperature is decreased below the KT transition, 
the amplitude of the flux noise, which is proportional to 
1/e Ril- l/e([jj = 0), rapidly decreases. The KT transi- 



V. CONCLUDING REMARKS 

In the present paper we have explored the fact that the 
typical experimental setup for measuring the magnetic 
flux noise for a superconductor, or JJA, corresponds to 
the case when the distance d to the pick-up coil is much 
larger than the relevant microscopic lengths. In this limit 
the flux-noise spectrum and the real part of the conduc- 
tivity are proportional. This proportionality seems first 
to have been anticipated in Ref. [?] and has also been ex- 
perimentally verified!!! In this paper we have studied this 
connection in some more detail. 

We have also demonstrated that both the shape and 
the characteristic frequency of the spectrum depend on 
the distance d to the pick-up coil. This means that no de- 
tailed conclusions can be drawn from simulations which 
presumes d = 0. Furthermore, there is a qualitative dif- 
ference between the vorticity-noise spectrum, which cor- 
responds to discrete events over a sharp boundary, and 
the magnetic- flux-noise spectrum which corresponds to 
spread-out objects over a boundary. The experimental 
situation corresponds to a spread-out magnetic flux and 
a pick-up coil at a large distance d, which is very differ- 
ent from some earlier simulations ^riwch calculated the 
vorticity-noise spectrum for d = 0.Bu Nevertheless the 
vorticity-noise spectrum for d = has a cj~ 3 ' 2 tail for 
higher frequencies which seems to match the experimen- 
tal results,Q'El whereas the magnetic-flux-noise spectra for 
d = does not have such a tail. In accordance with the 
present simulations, we suggest that the resolution of this 
dichotomy is that in the large-d limit the magnetic-flux- 
noise spectrum does have an intermediate region with a 
w -3/2 bghavjoj- anc j that it is this intermediate region 
which is seen in the experiments. 

The proportionality between the magnetic-flux-noisc 
spectrum and the real part of the conductivity implies 
that the noise spectra for a sequence of temperatures just 
above the KT transition should in a log-log plot have a 
common tangent with the slope —1. The existence of 
such ap-common tangent has also be verified in experi- 
ments,!^ as well as in the present simulations. We also 
explicitly demonstrated through our simulations that for 



small d there is no such common tangent. 

The existence of a common tangent is by itself not nec- 
essarily conclusive. For example, the experimental data 
for the JJA's in Ref. g correspond to the large-d case 
and the data have indeed a common tangent with slope 
— 1. However, the spectra at a fixed temperature seem 
to have a 1/oj behavior over a very large region, which 
differ markedly from the spectra obtained in our simula- 
tions for the RSJ model. 

The present simulations also suggest that immediately 
above the KT transition there should be a very small 
temperature region where the noise spectpun has an in- 
termediate interval with a l/u> behavior.E^I It has been 
suggested that the data in Ref. || might perhaps be re- 
lated— fo this temperature region closest to the transi- 
tion.^ However, at the moment there seem to be no acu 
cepted explanation for the 1/w-behavior found in Ref. g.0 

Finally, we showed that the amplitude of the flux-noise 
spectrum drops dramatically as the temperature is de- 
creased below the KT transition, and that at the same 
time the characteristic frequency increases. It should also 
be possible to observe this effect in experiments. 



ACKNOWLEDGMENTS 

The authors are grateful to P. Svedlindh and O. Festin 
for discussions of their experimental data and to D. Bor- 
mann for discussions of the theory. This work was sup- 
ported by the Swedish Natural Research Council through 
contract FU 04040-332. 



iu K.-H. Wagenblast and R. Fazio, JETP Lett. 68, 312 (1998). 

11 C. Timm, Phys. Rev. B 55, 3241 (1997). 

12 M. Capezzali, Ph.D. thesis, Universite de Neuchatel, 1998. 

13 A. Jonsson and P. Minnhagen, Physica C 277, 161 (1997). 

14 A. Jonsson and P. Minnhagen, Phys. Rev. B 55, 9035 
(1997). 

15 B.J. Kim, P. Minnhagen, and P. Olsson, Phys. Rev. B 59, 
11 506 (1999). 

16 P. Olsson, Phys. Rev. B 46, 14 598 (1992). 

17 Note that this only gives a measure of the flux. The mag- 
netic field associated with this flux is spread out; it is 
strongest in the vortex core region and decays on the length 
scale A. 

18 Note that the circulation of the normal current around a 
plaquette vanishes since it is proportional to ^ (9i —0j) = 
0. 

19 J.D. Jackson, Classical Electrodynamics, second edition 
(John Wiley & Sons, New York, 1975). 

20 M. Lax and P. Mengert, J. Phys. Chem. Solids 14, 248 
(1960); R. Voss and J. Clarke, Phys. Rev. B 13 556 (1976). 

21 In this respect our results and the ones in Ref. M agree com- 
pletely but differ slightly from Ref. M where an exponent 
somewhat smaller than 3/2 was obtained for the same case 
(i.e., w 1.2 at T > 1.1). 

22 In Fig. H there appears to be a slight upturn of the spectra 
at the largest frequencies. As in Ref. H, this is an artifact 
of the fast Fourier transformation used in the calculations. 

23 P. Minnhagen, in NATO ASI on Fluctuation Phenomena in 
High Temperature Superconductors, edited by M. Ausloos 
and A. Varlamov (Kluwer, Netherlands, 1997), p. 279. 



1 For a review see, e.g., P. Minnhagen, Rev. Mod. Phys. 59, 
1001 (1987); Proceedings of the 2nd CTP workshop on Sta- 
tistical Physics: KT Transition and Superconducting Ar- 
rays, edited by D. Kim, J.S. Chung, and M.Y. Choi (Min- 
Eum Sa, Seoul, 1993); Physica B 222 253-406 (1996). 

2 J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 
(1973); V.L. Berezinskii, Zh. Eksp. Teor. Fiz. 61, 1144 
(1971) [Sov. Phys. JETP 34, 610 (1972)] 

3 M.J. Ferrari, F.C. Wellstood, J.J. Kingston, and J. Clarke, 
Phys. Rev. Lett. 67, 1346 (1991). 

4 C.T. Rogers, K.E. Myers, J.N. Eckstein, and I. Bozovic, 
Phys. Rev. Lett. 69, 160 (1992). 

5 T.J. Shaw, M.J. Ferrari, L.L. Sohn, D.-H. Lee, M. Tin- 
kham, and J. Clarke, Phys. Rev. Lett. 76, 2551 (1996). 

6 P. Svedlindh and O. Festin (unpublished); M. Bjornander, 
J. Magnusson, P. Svedlindh, P. Nordblad, D.P. Norton, and 
F. Wellhofer, Physica C 272, 326 (1996). 

7 J. Houlrik, A. Jonsson, and P. Minnhagen, Phys. Rev. B 
50, 3953 (1994). 

8 I.-J. Hwang and D. Stroud, Phys. Rev. B 57, 6036 (1998). 

9 P.H.E. Tiesinga, T.J. Hagenaars, J.E. van Himbergen, and 
J.V. Jose, Phys. Rev. Lett. 78, 519 (1997). 





s 



(b) 



(c) 



FIG. 1. (a) The vortex dipole moment associated with the 
dual lattice point r is estimated by the circulating current 
around the plaquette. The magnetic field at the observation 
point x is the summation of the contributions from all such 
vortex dipole moments, (b) The I x I pick-up coil is separated 
by the distance d from the L x L square array of Joseph- 
son junctions, (c) The whole array is divided into quadratic 
enclosures. The plaquettes surrounding such an enclosure is 
denoted by S n (shaded area). The magnetic flux due the 
plaquettes belonging to S n is calculated. 
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FIG. 2. Comparison between the flux noise from the full 

calculation jEq. (|1J])J and from the approximate scheme 

[Eqs. (|l6|), (|lq), and (\L9j)]. The figure shows the flux noise as 
a function of time t, S(t, d) = ($(i, d)$(0, d)) for a 16 x 16 coil 
with the distance d — 5 from a 32 x 32 array at the tempera- 
ture T — 1.10 (in units of J/fcs). As seen the approximation 
scheme reproduces the full calculation very accurately. 
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FIG. 3. The weight factor q n (in arbitrary unit) for S„ plot- 
ted against the size 2n— 1 (in numbers of plaquettes per side). 
The data are for a 64 x 64 array with the coil size 32 x 32 
for different distances d. When the distance d between the 
coil and the array is very small, only plaquettes inside of the 
pick-up coil contribute to the magnetic flux, whereas all pla- 
quettes contribute for larger d. 
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FIG. 4. Comparison between the vorticity-noise spectrum 
Si and the magnetic-flux-noise spectrum S2 for d — 0. The 
data are for a 64 x 64 array at T — 1.1 with a 32 x 32 coil 
size. The vorticity spectrum has a <-o~ 3 X2 tail, whereas the 
magnetic-flux spectrum is closer to w -2 .Lj 
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FIG. 5. (a) The dependence of the flux-noise spectrum 
S(u>, d) on the distance d. The full drawn uppermost curve 
is the imaginary part of the dielectric function |Im[l/e(w)]|. 
The rest of the curves correspond to d = 0.1, 5, 10, 20, and 40 
(from bottom to top) plotted as u>S(u>, d) and the curves are 
shifted in the vertical direction for better comparison. The 
data are for a 64 x 64 array with a 32 x 32 pick-up coil at 
T = 1.1, except for d — 40 where a 128 x 128 array was 
necessary because plaquettes further away from the center 
contribute in this case. As d is increased u>S(u>, d) approaches 
|Im[l/e(oj)]|. (b) The frequency at the maxima for the curves 
in (a) are plotted versus the distance d. As d is increased, this 
frequency decreases and approaches the value for |Im(l/e(a>)|. 
(The full line is a guide to the eye.) 



FIG. 6. (a) The real part of the conductivity Re[<r(oj)] 
and (b) the flux-noise spectrum S(u>, d — 20) at tempera- 
tures above the KT transition. The curves for Re[cr(aj)] and 
S(uj, d = 20) have the same shape; for small frequencies they 
have a very weak oj-dependence, for somewhat larger u> there 
is an approximate uj^ 3 ' 2 -behavior, whereas for even larger uj 
the behavior approaches ui~ 2 . The curves in the log- log plot 
have a common tangent with the slope —0.9. 
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FIG. 7. Flux-noise spectrum in the small-d limit at tem- 
peratures above the KT transition [same as Fig. 0(b) but 
with d — 0.1 instead of d — 20]. In the case of small d we 
find neither a common tangent nor an appreciable range of 
w _1,5 -behavior. 
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FIG. 8. (a) The imaginary part of the dielectric function, 
|Imjl/e(u>)]|, and (b) the flux- noise spectrum multiplied by 
the frequency, uiS(u),d = 20), at temperatures above and be- 
low the KT transition. This again illustrates that both quan- 
tities behave in the same way. As the temperature is increased 
far above the KT transition, the maximum of [Im[l/e(o;)]| ap- 
proaches the limit value 1/2 (this corresponds to the Drude 
limit with e = 1) as denoted by the horizontal line in (a). The 
curves seem to develop a plateau as the KT transition is ap- 
proached from above, as is suggested by the T = 0.95-curves. 
As the temperature drops below the KT transition, the am- 
plitude of the flux noise rapidly decreases whereas the char- 
acteristic frequency increases, as is illustrated by the curves 
at T = 0.85. 
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